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Abstract 

We recalculate the next-to-leading order Altarelli-Parisi kernel using a method which relates 
it to the splitting amplitudes describing the collinear factorization properties of scattering ampli- 
tudes. The method breaks up the calculation of the kernel into individual pieces which have an 
independent physical interpretation. 
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1. Introduction 



The last quarter-century of experimental studies at colliders, complemented by theoretical investiga- 
tions, have taught us that perturbative quantum chromodynamics (QCD) gives an excellent description 
of the strong interaction probed at short distances. Indeed, the theory has reached a sufficient maturity 
that it is no longer the target of experimental studies, but rather a tool in the search for new physics 
beyond the standard model. 

Upcoming collider experiments at the Tevatron and the LHC will require continuing refinements and 
progress in our ability to make precise predictions of QCD and QCD-associated processes, along with 
an ability to give credible estimates of the associated uncertainties. Recent years have seen great 
progress in computing two-loop amplitudes |in|2l|3l in QCD, which is one of the ingredients essential 
to a next-to-next-to-leading order (NNLO) description of collider processes. 

The calculational framework for collider processes is based on our ability to compute short-distance 
matrix elements to increasing perturbative order in nonabelian gauge theories. It complements these 
process-dependent matrix elements with a general understanding of factorization, which separates out 
the process-independent long-distance aspects. The long-distance aspects of scattering processes are 
captured in the parton distribution functions of the scattering nucleon(s), and in fragmentation functions 
for identified outgoing hadrons. Up to subleading power corrections, such functions along with the 
strong coupling a, are the only ingredients needed from outside perturbation theory for a description 
of collider scattering processes. 

The parton distribution and fragmentation functions are functions of a momentum fraction, and of the 
scale Q^- at which the hadron is probed. As is well-known, their values at different 2^ are not inde- 
pendent, but are governed by the Altarelli-Parisi equation, whose kernel is computable in perturbation 
theory |4l|5l|6l0|8l|9l- It allows these functions to be evolved from a fixed scale to other values of 
(Jp-. The evolution kernels are known up to second order fl0l|Il]|l2l[l3[l4l[l5l[l6l|l7l[^T^. Curci, 
Furmanski, and Petronzio [10. IIJ used the light-cone gauge, whereas Floratos et al. obtained their 
results fT2lfT3lfT4lfT5lfT6ll using a covariant gauge. Subtleties related to the proper renormalization of 
the gluon operator in covariant gauges were later clarified by Hamberg and van Neerven l20l . Based 
on this work the second-order computation was rechecked more recently by Mertig and van Neerven, 
who also computed the polarized kernels at this order |21 1. The calculation of the polarized kernels 
was checked by Vogelsang [22 1 using the light-cone gauge. The use of light-cone gauge was also in- 
vestigated in refs. 03 1211 123- Its use beyond next-to-leading order is subtle |26| and not fully tested. 
An NNLO calculation is being undertaken, using the operator approach, by Moch, Vermaseren and 
Vogt 127J. Furthermore fits to the first moments of the evolution kernels [28 . .29 ] have already been 
used to determine the NNLO parton distribution functions l30ll3T]| . We believe it is useful to develop 
other methods applicable to an NNLO calculation, and that is our purpose here. 

Accordingly, we wish here to introduce an alternative approach to computing the kernel 1^ . We 
will do so using gauge-independent quantities describing the coUinear behavior of various amplitudes. 
This has the advantage of breaking down the computation into pieces which themselves already have 
a meaningful interpretation. It also allows us to avoid complications associated with prescriptions for 
light-cone gauge. A related analysis has been presented in refs. l33ll34l . In this paper, we focus our 
attention on flavor-independent contributions to the time-like gluon kernel. 

At present there exist two different methods to calculate the Altarelli-Parisi kernels. To contrast pre- 
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vious methods with the approach presented in the present work, let us consider the basic ingredients 
of factorization. The key feature of factorization in the QCD-improved parton model is the follow- 
ing: a hadronic cross section (more precisely, the leading-twist contribution) can be factorized into 
a hard scattering coefficient and parton distribution or fragmentation functions (so-called PDFs). All 
information about the hard scattering process is contained in the hard-scattering coefficient which is in 
the realm of short-distance physics and hence can be calculated entirely in perturbative QCD. Long- 
distance effects, including the matching of partonic states to hadronic ones and vice versa, are captured 
by the PDFs. These functions can be expressed as expectation values of composite operators between 
hadronic states. The matrix elements of the composite operators cannot be calculated in perturbative 
QCD, although they are (in principle) amenable to non-perturbative techniques such as lattice QCD. 
On the other hand, these non-perturbative matrix elements are universal, that is process-independent, 
and hence can be measured in one process and then used to obtain predictions for all other processes 
of interest. While the matrix elements themselves cannot be calculated perturbatively, their scaling 
behavior can be, and it is this scaling behavior which is captured by the Altarelli-Parisi equation. 
Their scaling behavior, or equivalently their anomalous dimensions, are determined by the ultraviolet 
singularities of QCD corrections to the matrix elements. 

The ultraviolet singularities, and hence the anomalous dimensions, are independent of the choice of 
external state in the composite-operator matrix element. For calculational purposes, we can therefore 
replace a hadronic state with a partonic one. This yields one method (the 'OPE approach') of calculat- 
ing the evolution kernels. 

The other method used in the literature is the so-called 'infrared approach'. Here one starts with 
unrenormalized quantities. The calculation yields singularities, in particular in the partonic, short- 
distance, cross section. (The singularities are typically regulated using a dimensional regulator.) One 
can distinguish between ultraviolet divergences and mass singularities, related to the coUinear emis- 
sion from initial-state partons (or from final-state partons fragmenting into identified hadrons). Soft 
singularities cancel between virtual and real corrections. Ultraviolet singularities are removed through 
the usual renormalization procedure, that is are absorbed into the definition of the physical coupling in 
terms of the 'bare' coupling. For mass singularities, the situation is more complicated. In the infrared 
approach, they are canceled by corresponding singularities in the 'bare' PDFs. Equivalently, they are 
absorbed into the definition of the physical PDFs in terms of the 'bare' PDFs. For this to be possible, 
the singularities must of course possess a universal, process-independent form. (They do.) 

The mass singularities in the PDFs determine their evolution with respect to changes in the reference 
scale, so the consistency of this approach means that the singularities determine the evolution kernels. 
In turn, a determination of the remaining singularities in the ultraviolet-subtracted hard-scattering co- 
efficient thus also yields a calculation of the evolution kernels. To determine these singularities, one 
could in principle calculate the complete partonic cross section for a specific process, and extract its 
singularities. In general, this is a formidable task. For example, in order to derive just the leading- 
order (LO) kernels one already needs to know the singularities of the next-to-leading order (NLO) hard 
scattering coefficient. 

Fortunately, this task can be simplified. One way is to use the observation of Ellis, Georgi, Machacek, 
Politzer, and Ross l35ll36l . that in axial gauge only a certain class of diagrams contribute to the sin- 
gularities. It is on this observation that the derivation of Curci, Furmanski, and Petronzio |10 11 1 is 
based. Their derivation is closely connected to the use of light-cone gauge whose treatment beyond 
next-to-leading order, as mentioned earlier, is not fully tested. 
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The method we shall use in the present paper is related but distinct. We use the knowledge of factoriza- 
tion at the amplitude level, more specifically the splitting amplitudes which describe the coUinear limits 
of scattering amplitudes. While it is intuitively clear that the factorization and universality of the 'mass' 
singularities in a hard-scattering cross section follow from this factorization at the amplitude level, up 
to now no expUcit demonstration has been given in the literature. We give such a demonstration, using 
the phase-space slicing method, and use the connection to derive the NLO evolution kernels. In the 
context of our derivation, the splitting amplitudes represent smaller parts of the calculation that can 
be verified independently. Furthermore, since they describe properties of on-shell amplitudes, they are 
manifestly gauge-invariant; the method we describe is accordingly independent of the use of light-cone 
gauge. 

Splitting amplitudes have been used extensively as a check on new calculations of amplitudes, since 
the constraint of obeying the correct behavior in all coUinear limits is a strong one. Indeed, the one- 
loop splitting amplitudes have also been used, via such constraints, to give conjectured forms 1.371 
for certain classes of one-loop amplitude with an arbitrary number of external legs (later proven by 
Mahlon EH), and as an aid in the initial derivation of another all-n class of amplitudes obtainable 
using the unitarity -based method l39l . 

In this work we will focus on the timelike g — > gg-kemel, with g denoting a gluon. To simplify the 
derivation we will consider a gauge theory with no light fermions, and with an additional massive 
colorless scalar coupled to the gluons via an effective vertex. 

At leading order, the splitting amplitude can be regarded as the amplitude for finding a parton of given 
momentum fraction inside a parent parton. It is therefore natural to expect the Altarelli-Parisi kernel, 
which describes the probability of finding a parton of given momentum fraction inside a parent parton, 
to be the square of the splitting amplitude. We make this correspondence more precise in section |3 
Beyond leading order, we expect that there will be virtual corrections to this picture. Indeed, they are 
given by the one-loop splitting amplitude, more precisely by its interference with the leading-order 
splitting amplitude. As usual in gauge theories, however, there are additional singularities due to the 
emission of soft or coUinear partons, and so we must also integrate over corresponding real-emission 
contributions. Heuristically, these two contributions, less the iterated leading-order kernel, give the 
next-to-leading order Altarelli-Parisi kernel. 

The outline of the paper is as follows. In section |2l we review briefly the factorization of color ordered 
amplitudes. In sectionQlwe setup the general framework. In the following section we illustrate the new 
approach by the rederivation of the leading-order kernel. We will discuss the leading-order derivation 
in great detail because parts of it will be reused in the derivation of the NLO kernel which we present 
in section |5] We give our conclusions in section |6l 

2. CoUinear Factorization 

The properties of gauge theories are easiest to discuss in the context of a color decomposition ll40ll4n 
02,,43j,44J. At tree level, for all-gluon amplitudes such a decomposition takes the form, 

4'\{ki,h,ai}) = £ Tr(r<'(i)..-r''("))Ar(a(l^',...,«^")), (2.1) 

aeS„/Z„ 
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where Sn/Zn is the group of non-cyclic permutations on n symbols, and j ' denotes the j-th momentum 
and helicity. We use the normalization 

Tr(r r*) = 5"^ (2.2) 

for the generators of SU(N). The same color decomposition as shown in eqn. (I2.lt holds for the 
amplitudes we shall consider, for the process (Sf ^ g---g, where (j) denotes a colorless heavy scalar. 
One can write analogous formula for amplitudes with quark-antiquark pairs. The color-ordered or 
partial amplitude A„ is gauge invariant. In the coUinear limit, ka \\ h, of two adjacent legs, the color- 




Figure 1 : A schematic depiction of the collinear factorization of tree-level ampli- 
tudes, with the amplitudes labeled clockwise. 



ordered amplitude A„ is singular. (It is finite when the two collinear legs are not adjacent arguments, 
that is when they are not color-connected.) This singular behavior has a universal form expressed by 
the tree-level factorization equation, 

Ai°)(l,...,.^^^^^...,.) 

g, £ Split':;^(a^^^^'')Al°J^(l,...,(a + Z7)^.. .,«) + •••, (2.3) 

ph. pol. o 

where the dots represent terms which are finite in the limit. In this equation, Split''^'^'^ is the usual 
tree splitting amplitude, and the notation 'a + b' means ka +ki,. The QCD couphng is denoted by 
gs- The notation 'ph. pol.' indicates a sum over physical polarizations only. ('Physical' here is in the 
sense of 'transverse', and their number may depend on the number of dimensions and on the variant 
of dimensional regularization employed.) This factorization is depicted schematically in fig.Q At tree 
level, one may derive the sphtting amplitudes from a string representation |45| or from the Berends- 
Giele recurrence relations ll46l . It is characteristic of gauge theories that the splitting amplitude has a 
square-root singularity. Split ~ 1 / ^Jsab, rather than a full inverse power of the two-particle invariant 
Sab- Similar formulae hold in the triple-coUinear case [47 48 , 4^ I5U1 . 

At one loop, the color decomposition analogous to eqn. ( 12.11 ) is 

L«/2J + 1 

An{{ki,h,ai}) = Y,nj £ £ Gr„;, (a) AiiUa), (2.4) 

/ C=l <5eS„/S,r,, 

where [xj is the largest integer less than or equal to x and is the number of particles of spin /. The 
leading color-structure factor, 

Gr„;i(l) =A^ Tr(r"' •••r"") , (2.5) 
is just N times the tree color factor, and the subleading color structures are given by 

Gr„.^(l) = Tr(r" • • • r'-') Tr (r'-- • • • T"") . (2.6) 
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Figure 2: A schematic depiction of the colUnear factorization of one-loop amplitudes. 



Sn is the set of all permutations of n objects, and S^-c is the subset leaving Gr,,;^ invariant. The de- 
composition eqn. (I2.4t holds separately for different spins circulating around the loop. The usual 
normalization conventions take each massless spin-7 particle to have two helicity states: gauge bosons, 
Weyl fermions, and complex scalars. (For internal particles in the fundamental (N +N) representation, 
only the single-trace color structure (c = 1) would be present, and the corresponding color factor would 
be smaller by a factor of N.) 

The subleading color amplitudes An-c>i are in fact not independent of the leading color amplitude 
Rather, they can be expressed as sums over permutations of the arguments of the latter 
I5T1 . (For amplitudes with external fermions, the basic objects are primitive amplitudes |52| rather 
than the leading color one, but a similar dependence of the subleading color amplitudes on the leading- 
color ones holds.) As a result, it suffices to examine the collinear limits of leading color amplitudes. 
The collinear limits of the subleading color then follow using this relation. The leading color one-loop 
amplitudes obey the following factorization 15111531 . 

£ Lsplit':;^(a^",Z7^'')Ai^2i(l,...,(a + Z7)^...,«) 

ph. pol. a \ 

+gl Spliti-^''(a"",^"*)Al"2i(l, . . . , {a + br, ...,«)). (2.7) 
This factorization is depicted schematically in fig. El 

This form was originally deduced from explicit calculations of higher-point amplitudes l5Tll52ll . but 
can also be proven more generally using the unitarity-based method l53l . The latter proof also provides 
an explicit formula for the one-loop splitting amplitude. We used it l54ll to calculate all the one-loop 
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splitting amplitudes relevant in QCD to all orders in the dimensional regularization parameter 



z = {4-d)/2 (2.8) 

with d being the dimensions of the spacetime. A subset of the terms higher order in e are needed for 
singular phase-space integrations in NNLO jet calculations |55|. Bern et al. L56. .57. . 58 1 derived the 
one-loop splitting amplitudes from an analysis of one-loop integrals. 



3. Framework 

In this section we derive a relation between the singularities in a 'partonic' (unsubtracted) cross section 
and the evolution kernels. We use this relation in following sections to compute the leading-order 
and the next-to-leading order kernels. Our derivation follows closely the one given in ref. |59|. As 
mentioned earlier we restrict attention to pure gluonic QCD together with an uncolored massive scalar 
(|) to which the gluons couple via a higher-dimension operator. Such a coupling could, for example, be 
induced via a heavy quark loop which is integrated out. A similar derivation would of course hold for 
a gauge theory with fermions as well. 

Consider the production of a glueball G in the decay of the massive scalar (|). In order to calculate 
the corresponding decay rate or differential distributions thereof, we need as ingredients the subtracted 
decay rate T^{(Sf ^ g+X) for the production of gluons in the decay of the massive scalar, and the 
gluon-to-glueball fragmentation function 'D^^q. Using these quantities the energy distribution of the 
daughter glueball is given by the following expression. 



dr{(^^G+x) _ dr'^{<^^g+x] 



dxc dxn ^ ^ 



where the convolution '(g)' is defined by 

[f^g]{x) = C f dydzf{y)g{z)h{x-yz). (3.2) 

JO JO 

The energy fractions xq and Xg are normalized to the mass of the scalar particle: 

XG,g = (3.3) 

where Eq [Eg) denotes the energy of the glueball (gluon). We use the subscript R to indicate that both 
the fragmentation function as well as the subtracted decay rate depend on the subtraction method used 
to define them; we have not put in an explicit argument to show the dependence on the factorization 
and the renormalization scale. The universality of factorization allows us to write down a formula for 
the 'bare' decay rate r^((|) ^ g + X) of the scalar (|) into an identified gluon in terms of the subtracted 
decay rate: 

.r>-.,+x)^.r>-.g+x)^^^^^ (3.4) 

The 'bare' decay rate is unphysical, as it describes the 'probability' of finding a gluon with a given 
energy fraction inside a jet, while of course a lone gluon is not a colorless physical state. It will thus 
contain coUinear or 'mass' singularities, as will (Dg^g. Our purpose in considering r^((j) ^ g + X) is 
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precisely to extract these singularities. In the following we will assume that all singularities (ultraviolet 
(UV), soft, and mass singularities) are regulated via a dimensional regulator. The advantage of the bare 
decay rate r^((t) ^ g + X) is that it can be calculated purely in perturbative QCD. As suggested by 
eqn. IB. 41 the mass singularities which are present in "Dg^g must exactly match those in the (unphysical) 
partonic decay width r^((j) g + X). To see this, expand 'Dg^g in a^. 



'Dg^giz) = 5(1 -z) + ^^^%iz) + {^y4%iz) + 0(a.^) 



and then invert eqn. (13 .4t to obtain. 



I® 



-1 



+ 



g^g ' 



I'D 



(1) 

g^g ' 



(2) ( 
g^g\ 



z)]+0(a' 



where we have used the identity 

5(1 -z) = 



5(l-.).|.i2,(z).(|).fi.(z) 



(2) 
g^g 



iz)) 



(3.5) 



(3.6) 



+ Oial) (3.7) 



to invert ©g^g. The left hand side of eqn. (13.61) is a finite quantity, and thus the right hand side must 
be so as well. If we now expand the 'bare' partonic decay width in a^, we see that order by order 
the singularities in the partonic decay width must be canceled by those which appear in "Dgl^g. In 
particular, using 



(3.8) 



we obtain 



dT^i^^^g + X) 



dx„ 



271 



1® 



g^g 



+ 



g^g^ 



1® 



(1) 

g^g " 



^^%{z] 



+ 0(a?). (3.9) 



(The reader may worry about powers of implicit in the higher-dimension operator coupling gluons 
to the heavy scalar. Such factors are frozen at the scale where the operator is generated, and in any 
event do not enter into the following arguments.) The usual application of this expansion is to the 
calculation of the differential decay width <ir^((|) —>■ g + X)/dxg: start with the 'bare' partonic decay 
width h^'^ and use eqn. (13. 9t (after ultraviolet subtractions as well) to obtain the finite subtracted decay 
width, which predicts the decay of the scalar (|) into a glueball or more generally into hadrons. Here, we 
will use our knowledge of the collinear divergences in the 'bare' partonic decay width r^((|) —>■ g+X), 
along with eqn. (I3.9t . to determine 'Dg^g. 
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Through o(a^), we can express Dg^g in terms of the evolution kernels P^'^^ and 

5(1 -z) 



271 £ 



z 



+ 



271^ ^ 



JLp(o)^p(o) + JLpp(o)_iip(i) 

2£2 4£2 2 £ 



+ 0{al). 



As usual Po denotes the first coefficient of the QCD (3-function, 

„ 1 , , Hf=0 11 

P()= 3(llA^-2«/) '= —N, 
and Sg is the usual factor appearing in the modified minimal subtraction (MS) scheme: 



(3.10) 



(3.11) 



(3.12) 



with Y the Euler constant. The general structure of this expansion follows from the renormalization- 
group equation in the MS scheme, which we choose. In particular, using eqn. (I3.10t we obtain for the 
fi dependence of 2)f_^G- 



-^g^g^[ms,e)^%ig) 



g^G 



g^G 



(3.13) 



where we have used the P-function in d dimensions P(av,£) in the MS scheme: 

p(a„£) =f?—a, = -a,(£ + ^5Epo) + o(a,^). 
Inserting eqn. (I3.10t in eqn. ( 13. 9t we obtain for the subtracted differential decay rate. 



(3.14) 



dx^ 271 



2£2 



4£2 



/j(2)+/,(l)^l5^p(0)(,) 

(3.15) 



2£ 



The subtracted decay rate r^((|) g+X) calculated via eqn. (I3.15t is the MS-subtracted decay rate. In 
the following two sections we illustrate how to calculate the divergences in h^^^ and h^'^\ and thereby 
determine the kernels P'"' and P^^\ 



4. The Leading-Order Kernel 

In order to compute the leading-order kernel, P^^\ we must isolate the collinear singularities in h^^\ In 
terms of the matrix elements, the leading-order partonic decay rate is given by the following formula. 



/i(0)(. 



I dR^hM) W''\'^^g{h)g{k2))\'mxg-x,) + h{xg-X2)] (4.1) 



8 



where dR^{k\ ,kn) denotes the phase-space measure in d dimensions (including the symmetry factor 
l/«!) for « gluons with momenta ki,...,kn, 

The leading-order amplitude for the decay (|) ^ gg is given by . The energy fraction x,- of the 
identified gluon / is defined as in eqn. (I3.3t but with Ec^g replaced by its energy £",. At next-to- 
leading order the partonic decay rate gets two additional contributions, from virtual and real-emission 
corrections, 

■g(^l)g(^2))^('"*((^-g(^l)g(^2))] 
■^2)], 

g{h)g{k2)g{k,))\^ 

X2)+6(Xg-X3)], (4.3) 

where A ((j) gg) is the one-loop amplitude for the decay (|) — > gg while ggg) is the ti^ee- 

level amplitude for the real-emission process (|) —>■ ggg. 

Both contributions in eqn. (14. 3t contain collinear singularities associated with the identified gluon and 
will thus contribute to the kernel /'("^(z). (The individual contributions also contain soft and other 
collinear singularities.) In the virtual corrections h[P the singularities arise from the loop integration. 
The general structure of infrared divergences of one-loop amplitudes is known, we may extract them 
without further knowledge of the specific process at hand. This contribution will only contribute to the 
term proportional to 8(1 — z) in P'"^ (z). 

Let us then turn to the contribution from real emission which will determine the 'non-trivial' z de- 
pendence of P^^\z)- (The 5(1 — z) contribution can always be determined from the z dependent part 
through appeal to various sum rules, for example that associated with momentum conservation.) In 
the real emission contribution the singularities arise from the phase-space integration of singular 
terms in the matrix elements. These may be extracted by the phase space-slicing method l60l l6n . 
We will not need the whole apparatus developed for cross-section calculations, but only the following 
basic elements. The basic idea in the phase-space slicing method is to split phase space into resolved 
and unresolved regions. In resolved regions, all outgoing partons are 'resolved', which is to say none 
become soft or collinear. The unresolved regions are the remaining regions of the phase space; in these 
regions one or more final-state partons may be soft, or one or more pairs may become collinear. 

For the three-gluon final state, we may use the following partition of unity to separate the various 
regions, 

1 = [0(512 -i'mm)+0('5'mm-'5'l2)] [0(^23 -■5'min)+0('5'min-'5'23)] 
[&{sii, - Smin) + &{Smm " S13)] 
= @{S12 - Smm)@{s23 " Sniin)@{si3 - ^min) (1, 2, 3 hard) 
+ ©(■Smin - ■512)0(^23 " ^mm)0(^13 " ^min) (1, 2 CoU.) 
+ 0(^12 - Smm)@{Smm " S23)@{si3 " ■^min) (2, 3 COll.) 



g/.(i)(x,) = '^{hi'\xg)+h9\xg)), 

—h[^\xg) = f dR''{kuk2)2Re[A^^\(^^ 
X [5{xg — xi ) + 5{xg — 

^hi'\xg) = I dR''{k,,k2M)\^'^°\<^^ 
X [6(x^ — xi ) + 5{xg — 
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+ 


0(512 - 




'^min)0('^min 




(1,3 coll.) 


+ 


©(■^min 


-512)0(523 - 


" '^min ) (■^min 




(1 soft) 


+ 


®('^min 




-523)0(513 - 


~ '^min) 


(2 soft) 


+ 


0(512- 


~ '^min)0('^min 


-523)0(5min 


-^13) 


(3 soft) 


+ 


®('^min 




-■^23)0('«rain 




('double unresolved') 



(4.4) 

with Sij = 2ki • kj. The arbitrary parameter 5rain introduced in eqn. (I4.4t controls the boundary between 
resolved and unresolved regions. In parentheses we have classified the different contributions into 
resolved, soft and collinear contributions. The 'double unresolved' contribution does not contribute 
because it is kinematically forbidden. Thus only the contributions classified as soft or collinear will 
yield singularities. 

As we shall see, in fact only the collinear contributions survive at the end after combining virtual and 
real-emission contributions. As one would expect, the singularities arising from soft regions cancel 
between these two types of contributions. The evolution kernels are thus determined solely by the 
collinear or 'mass' singularities. From a practical point of view, such cancellations are an important 
consistency check on the calculation. 

The cancellation of soft singularities can be seen heuristically from the form of the partonic observable 
we are 'measuring'. For an n-parton final state, it takes the form, 

n 

On = ^5(Xg-.X,). 

i=l 

We see that 0„ — > 0„_i when one gluon becomes soft, but this does not happen when a pair becomes 
collinear. That is, the observable is soft-finite but not free of collinear singularities. Were the 5- 
functions not present in our integral, the sum of the virtual and real would of course be finite because 
it would just be the NLO correction to the total decay rate. With the 5-functions present, but one of 
the gluons soft, the corresponding term does not contribute, and all other terms are insensitive to the 
soft momentum. We can thus integrate the real-emission contribution over the soft region, generating 
poles in £ that cancels the corresponding singularity in the virtual corrections. In the case of collinear 
gluons, however, the two 5-functions depending on the momentum fractions of the collinear gluons 
will prevent us from integrating over the appropriate region of phase space. The collinear divergence 
in the virtual corrections will therefore not be fully canceled. The left-over singularity is precisely the 
term we are trying to extract. 

We can use the symmetry of the matrix elements and of the phase-space measure to restrict attention 
to the contribution where gluons 1 and 2 are collinear. The other two collinear contributions follow 
through symmetrization. Using the color decomposition introduced in section |2l we may write, 

^^'\<^^gig2g3)= £ Tr(r<'(')r<'P)r«<'ffl)A«')(a(l),a(2),a(3)). (4.5) 

aeS3/Z3 

Calculating the square of (")(([) g 18283) yields different color structures. For the determination of 
P^^\z), it is sufficient to consider the leading-color structure (denoted by the subscript 'Ic'): 

\^^'H'^^8i8283)\\ = Tr(rr^nTr(rr^r) £ \A^'Hi,j,k)\' 



(01)=(123),(132) 



Ic 



A^3(|a(o)(1,2,3)P + |a(")(1,3,2)|2). (4.6) 
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Within these terms, examine the contribution from the phase-space region where gluons 1 and 2 become 
collinear, 

J dF^ikxMM) \^^^\^ g\g2g?.)\\ [6(xg -xi) + 8(xg -X2) + 5(xg -X3)] 

X [5{Xg - ZXp) + 6(Xg - ( 1 - Z)xp) + 6(Xg - X3)]0(i'min - Su)®{s21, " Si^m)®{sn " ■^min) • (4.7) 

The factor 2 in front accounts for the two different color orderings. In deriving eqn. (I4.7t we have used 
the factorization of the amplitudes as discussed in section |2l Using these results the factorization for 
the squared matrix element is given by 

where we have eliminated cross terms of the form Split'^*^^ {lij^' x (^Split^'^'^ {kf' ,k^^)j . These terms 
vanish upon azimuthal integration. We have also made use of the factorization of the phase-space 
measure in the collinear limit L60J . 



dR'ihMM) ^ ^-dR\PM)dRt,,,{k,M) (4.9) 

where 

dRiiiihM) = ^{sminfdsndz[snz{l-z)r'' (4.10) 

and 

' 1 (15^)', (4.11) 

167l2r(l-£) V^min / 

We have defined z to be the momentum fraction in the collinear limit, 

ki=z{ki+k2)=zP, and k2 = {I - z){ki +k2) = {I - z)P (4.12) 
We work throughout in = 4 — 2s dimensions. The leading-order splitting amplitude is given by 

\/2 

Split'''='=(l,2) = (-£i -sj/ci-ep +/c2-eiep-e2 - ki ■ £2^1 ■ £p) ■ (4.13) 

If we are only interested in the unpolarized kernels we have to sum over the final gluon polarizations 
and average over the fused-leg polarizations, 

£|Split-'=(l,2)p^i £ |Split;^-(l^>,2^^)|2 = Al£!_i±I)!^Ap(,). (4.14) 

The averaging depends in general on the variant of dimensional regularization, this form holding for 
the conventional scheme Ii62il . Inserting eqn. (14. lOt and eqn. (14. 141) in eqn. (I4.7t we obtain 

1,(1 -,),-.|.(0)(P,,,)|V(,) 

3e J Jz{h,p) 

X [5(Xg - ZXp) + h{Xg - (1 - Z)xp) + 6(Xg - X3)]0(53p - imin), (4. 15) 
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with 

f(/,j) = — • (4.16) 
^'7 



eqn. (I4.15I) has almost the form we want: if we extend the region of the z-integral by adding and 
subtracting 

[,R.:iP.,,)(r\ f )</d.(i-.)i-i^'°'(/>.*,)ivw 

X [h{Xg - ZXp) + 5{Xg - (1 - Z)xp)]&{s3p - 5min), (4.17) 

we can write the expression in eqn. (14. 15t as a convolution plus an additional term: 
2 



where 



+ -A^l ^//?f(l,2)|^(°)(l,2)|2[8(x,-xi) + 6(x,-X2)]c('"(l,2,l), (4.18) 
= \6K^9i,N{8{l-zH-2{z{l-z)y^p{z)), (4.19) 



^ = C dzzr\l-z)-^p{z) 
Jo 



3(l-£)(4-3£) r^(l-£) 
'2£(3-28)(l-28) r(l-28) 



2 11 1 , 67, , , 



and 



C^'HiJ,k) = -2^ r'^''^dz[z{l-zT'p{z). (4.21) 



The other coUinear regions - in which || ^3 or /ci || - yield the same result. In the sum of the three 
we get thus a factor of 3 which cancels the factor 1/3 from the phase space measure. 



Using the factorization of color-ordered amplitudes in the soft limit L45.I 

A(°)(l,...,/^^^,/^...,«)■^'gsA^soft'-'=(/,^,70xA(o)(l,...,/^^/^...,^) (4.22) 

together with the factorization of the phase space measure in the soft limit 1601 . 

dR''{i,sJ) ^dRU{i,s,j)dR''{i,j) (4.23) 



with 

dRU{i,sJ) = 9{,{s^^r ('-^) " (4.24) 

V ^'7 / ^'7 

we obtain (after relabeling) 

X0(5min -■512)0(^23 " ^min)0(^min " ■^is) + (2 SOft) + (3 SOft) 

= y"j/?^(l,2)[5(x^-xi) + 6(xg-X2)]|^<'''(l,2)|22A^5(°'(l,2) (4.25) 
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for the singular contribution from the soft regions. The soft factor 5 (1 , 2) is given by 

5(")(1,2) / ^/<ft(2,3,l)g2|Soft''--(2,3\l)|20(s„i„-^23)0(^min-^i3). (4-26) 

As we shall see, we will not need the explicit result for Soft'^'''^. Combining the soft and coUinear 
contributions we obtain 



Ic, sing. 



a.s 1 

271 £ 



+ J J/?''(l,2)[5(x,-;ci) + 5(;c,-;c2)]|;i(°)(l,2)|22A^(5(")(l,2) + c(°)(l,2,l)](4.27) 

for the singular part of the real emission contribution. To complete our analysis of the singular part in 
/i^^) it remains only to add in the contribution from the virtual corrections. The color decomposition of 
the one-loop amplitudes in pure gluonic QCD is given by (c.f. section|2ll 



1 •••§«)= L Gr„;c(a)A„;f(a(l),...,a(«)), 

ceS„/S„-c 



(4.28) 



with Gr„;c(l) defined in eqn. (12. 6t . Once again we are only interested in the leading color-structure. So 
it is sufficient to include Gr„;i (l)A„;i (1, . . . ,«) = N Jr (T''' ■ ■ ■T"")a[1\\, . . . ,n) in our analysis. The 
singularity structure of the one-loop color-ordered amplitude a|J^(1,2) is known l60l l63l . Alterna- 
tively, we may reason as follows: suppose we are calculating the real-emission contribution to the total 
(|) decay rate. In this case we must replace the sum over 8-functions by 1 in the derivation above. We 
see immediately that the singular contribution of the soft and collinear regions to the total decay rate 
reduces to 

' J/?^(l,2)|.^ (")(!, 2)|22Af(5(°)(l, 2) + ^("'(1, 2,1)). (4.29) 



On the other hand, unitarity 1641 1651 dictates that the total decay rate must be free of soft and collinear 
singularities, order by order in perturbation theory. (Note that the ^min -dependent terms on the right- 
hand side are finite.) This implies that the singular contribution from the leading-color one-loop am- 
plitude satisfies 



^2Re(AfJ)(l,2)A('')*(l,2)) 



sing. 



-|A(0)(l,2)|22A^(5(°ni,2) + c(°)(l,2,l)) 



(4.30) 



sing. 



As a consequence we see that the singular contribution from the virtual correction cancels the second 
term in eqn. (I4.18t . We thus obtain 



Ic, sing. 



(4.31) 



for the surviving singular term. Note that h^'^^ should not be expanded in £. Comparing with eqn. (13.151) 
we find 



p(o)(,) = -ir«')(z) 
where we have used 



= 2n(^5{1-z) + - + 

d=4 V 12 Z 



1 



1-z 



(1-2) 



-1-E 



1 



5(1 -z) + 



1 




"In(l-z)" 




-£ 

+ 


1-z 



■z^+z-2 



+ 0(£2). 



(4.32) 



(4.33) 
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In these equations, the plus prescription defines distributions via, 

[F(z)]+ = lim|0(l-z-ri)F(z)-8(l-z-ri) '^F{y)dyY (4.34) 
so that if g{z) is well-behaved at z = 1, then 

/'*(r^ = jf'«M^+s(,)l„(,-.-). (4.35) 
■*,„i[14i^V = /' feW-V'"°"-'-> +^ln^(l-.-). (4.36) 



1-z 



1-z 2 



The phase-space slicing used here shares many features with that used for the calculation of a typical 
infrared-finite observable. In contrast to those observables, which necessarily allow the recombination 
of coUineai" partons in the real-emission terms, and hence allow integrating over coUinear phase space 
in a process-independent way |60|, the partonic decay rate r^{<^ g+X) describes an unphysical 
'probability' of finding a gluon with a given energy fraction inside the jet. This does not allow the 
recombination of two collinear gluons to a hard one, and a further convolution with a physical state 
distribution function is necessary to produce an infrared-finite observable. The uncanceled singularity 
does have a universal form, however, which is why this 'unphysical' object can be used to extract 
it. The singularities associated with soft-gluon emission do cancel, because too soft a gluon won't 
contribute to the energy fraction of any final-state hadron, and hence will drop out of the calculation. 



5. Next-to-Leading Order Kernel 

In order to calculate the NLO kernel we must study the singularities in the NNLO decay rate. Again 
keeping only the contribution dominant in the number of colors (denoted by the subscript Ic) the decay 
rate is, 



NNLO, Ic 



N' f dR'{l,2) 2Re(A;,')(l,2)A(«>(l,2))£5(^,-^,) 

+ N' f dR\\,2) |A;j)(l,2)|2f 5(x,-x,-) 

i=i 

+ 2N' I dR\\,2,3) 2Re(A(p(l,2,3)A(")*(l,2,3))£6(x,-x,-) 

+ 6N^ / J/?^(l, 2,3,4) |A(o)(l,2,3,4)|2£8(xg-xO. (5.1) 

i=i 

The factors of 2 and 6 in front of the integrals in the last two terms are combinatorial. Once again we 
are interested only in the singular terms. The singularities in the first two terms will contribute only to 
the coefficient of 5(1 —z) in the NLO kernel. The structure of the singularities, needed for the direct 
computation of this coefficient is now known in a general two-loop amplitude |66,,67,,68J. However, 
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the coefficient of the 8-function can always be computed using the sum rule constraints on the kernels. 
We will therefore not compute it directly, and turn our attention to the remaining terms in the kernel. 
These are determined by the three- and four-parton final states, which we discuss in the following two 
subsections. 



5.1. Three-Parton Final State 

The treatment of this contribution is very similar to the treatment of the three-parton final state in the 
computation of the leading-order kernel. In particular, we may use the same phase-space slicing. Here, 
we must consider the collinear limits of one-loop amplitudes in addition to those of tree amplitudes. 
Thus in addition to Split'"^'^, the one-loop splitting amplitude Split^"'°°P also makes an appearance. 
Symmetry again allows us to focus on the configuration k\ || ^2; the other two collinear regions will 
give equal contributions. We are interested in the singular contribution arising from 

2 / J/?''(l,2,3) £6(^,-x,)2Re(A}J)(l,2,3)A(o>(l,2,3)) 
1=1 

X0(^min - Sn)&{s23 " Smm)@{s\3 " ■^min)- (5.2) 

Using the factorization of the phase-space measure eqn. (I4.9t . along with the factorization of the lead- 
ing color one-loop amplitudes iBTIB^ (cf. eqn. ilJl ). 

A[l\\,...,i^-J'~^,...,n) ^ g,Split'^{i^-j'~n>^Al^\l,...,ii + j)\...,n) 

+gl Splitit °P (/^' , /O X a(0) ( 1 ,...,(/ + ;)",...,«) , (5.3) 

we obtain 

/ dR\p,3) / dzw-z)rpm{^,-zzp)+K^,-{\-z)zp)] 

3e J Jz{3,p) 
x2Re[A[l\p,3)A^-°>{P,3)] 

— ^^min-^ J dR'iP,3)J_^^^^ dz [z{l -z)]-'[5{x,-zzp) + 5{x, - {l-z)zp)] 
X 2|a(°) {P, 3) p£Re[(^i2) Split'"'™? (1 , 2) Split'^'='=* (1 , 2)] 

pol. 

+^ J dR\P,3) 5(x,-X3)c(°)(3,P,3)2Re(AjJ'(P,3)A(°)(P,3)*) 

+^|j/?^(P,3)5(x,-X3)c(')(3,P,3)|a(0)(P,3)|2, (5.4) 



where 



C^'\iJ,k) = -^s^,-^ r'^^^''\z[z{l -z)]-'=2£Re[(.i2)'+^Split'-'°°P(l,2)Split'--(l,2)]. 

(5.5) 

Note that ^ 

^Re[(5i2) '+*= Split'-'™P (1 , 2) Spht'^'^'^* (1,2)] (5.6) 

pol. 
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does not depend on su- Extending the region of the z-integral (by adding and subtracting the corre- 
sponding term), neglecting 5(1 — z)-type contributions, adding the contribution from k2 \\ ^3 and k\ \\ kj, 
and relabeling, we obtain 

^N^J dz J dR'{\,2) [5(x,-zxi) + 5(x,-zx2)]2Re[A(i)(l,2)A(0)*(l,2)]l/:(0)(z) 

i^y^^I 1 ^^''(l'^) [6(x,-z;ci) + 5(xg-zx2)]|A(°)(l,2)|2l^(^)(z) 



+ 



where 



pol. 



The explicit results for Split'"'°°P(l,2) can be found in refs. I5ni57ll58ll54l: 

Spliti-^°°P(1,2) = n(z)Split"--(l,2) + r2(z)^^^|^ (512S1 -82 -2^2 -ei^i -£2) 



with 



ri(z) 



2 V -^12 



[z/i(z) + (l-z)/i(l-z)-2/2], 

/2, 



(1- 2s) (3 -2s) 



(5.7) 



-(167l2)2^.^f„A^2(^(j_^))-e£R^((^^^)l+e3plijl-loop(i^2)Split'^-(l,2)*). (5.8) 



(5.9) 



(5.10) 



fiiz) 



-r(i -s)r(i +s)z- 1-^(1 -z)^ - - + -^^ 



-r(l - e)r(l +s)z-i-^(l -z)^ - - + ^-^—^ + -Li2(z) 

z z z 



2Fi(s,s;l+s;z) 

+ 0(s), (5.11) 
(5.12) 



where 2^1 is the Gauss hypergeometric function, Li2 the dilogarithm, and the standard one-loop pref- 
actor is, 

r(i+s)r2(i-s) 



(5.13) 



Using the above results we find that 



n (z) = Hjc 



1 1 



1 



1 



-^12/ V £ £ 

and thus for the virtual contributions to K^^\ 



+ - (ln(z) + ln(l - z) ) - - ln( 1 - z)" + ln(z) ln( 1 - z) - - ln(z 



(5.14) 



K^'iz) 



1 2 



{16ti^9{,YN^{ p{z) + - (ln(z) + ln(l -z 



-21n(l - zf - 21n(z) ln(l - z) - 21n(z)2 + ^Ti^) + \\ + 0(s) . (5. 15) 
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The above equations give the unrenormaUzed spUtting ampUtude; the renormahzed one is, as given by 
eqn. ( 15^ - eqn. (15m 



which adjusts eqn. ( I5.15t by, 



3 Ve 

,2 N -e 



167C2 £ 6 



ln(z)-In(l-z) U(z) + C>(£) 



(5.16) 



11 /I 



3 Vs 



In(z)-ln(l-z)U(z) + (9(£). (5.17) 



5.2. Four-Parton Final State 



We turn next to the four-gluon final state. The slicing of phase space is now more complicated than in 
the three-gluon case. A new feature arises: we can have double unresolved contributions. Such contri- 
butions can originate from two soft gluons, one soft gluon and a collinear pair, from two independent 
coUinear pairs or from a triple collinear configuration. To derive a suitable slicing we start with the 
following partition of unity, 

1 = [0('yi2-5min) + 0(^min-^12)][0('y23 -5min)+0('Smin-^23)] 

X [0(.S34-'5rain) + 0('5rain-'534)][0('5l4-'5min)+0('5min-'5l4)]- (5.18) 

Each term in eqn. (15.11) is expressed in terms of a single color-ordered amplitude, and hence the van- 
ishing of an invariant Sjj for non-adjacent gluons /, j does not yield a singular contribution. Expanding 
the right hand side of eqn. (I5.18t we get a decomposition of the four gluon phase space into sixteen 
different regions. We can classify these as chown in tabled 

In the first type of region, all invariants would be smaller than Sy^^^, but this is kinematically forbidden, 
and so it gives rise to no contribution. In the last type of region, all invariants are greater than ^min. 
so that no singularities arise; we can set aside this region as well. In the penultimate type listed, one 
nearest-neighbor invariant is smaller than s^i^ while the other three are greater than it. These give 
rise to single-unresolved contributions for a collinear pair of gluons. These regions give the NLO 
corrections to processes with three distinct jets. Their treatment follows that of the single-unresolved 
regions in the three-parton final-state in section ID As given, they contribute only to the leading-order 
kernel (and could be used to compute it had we not already used the NLO decay rate to do so). 

It will nonetheless be convenient to add and subtract terms as done in section |4] for the three-parton 
final state; the subtracted terms will then enter into the calculation of the NLO kernel. (This amounts 
of course to shifting terms from other regions.) The contribution from the 1 || 2, for example, is given 
by 

r ^ 
6A^M J/?^(l,2,3,4) |A(°)(l,2,3,4)|2£5(xg-x;) 

= f /j/?^(l,2,3) /^'^z|A(«)(l,2,3)p|V)(z) 
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Sl2 


S23 


S34 


J4I 


Comments 


< 


< 


< 


< 


Triple-unresolved: kinematically forbidden. 


< 


< 


< 


> 


1, 4 hard & separated. 1 | 2, 3 soft; or 3 | 4, 2 soft; or 2, 3 soft. 


< 


< 


> 


< 


3, 4 hard & separated. 2 | 3, 1 soft; or 1 | 4, 2 soft; or 1, 2 soft. 


< 


> 


< 


< 


2, 3 hard & separated. 1 | 2, 4 soft; or 3 | 4, 1 soft; or 1, 4 soft. 


> 


< 


< 


< 


1, 2 hard & separated. 2 | 3, 4 soft; or 1 | 4, 3 soft; or 3, 4 soft. 


< 


< 


> 


> 


1, 3, 4 hard. 1 || 2 || 3; or 1 || 3, 2 soft; or 1 ^3, 2 soft. 


< 


> 


> 


< 


2, 3, 4 hard. 2 || 3 || 4; or 2 || 4, 3 soft; or 2 ^4, 1 soft. 


> 


> 


< 


< 


1, 2, 3 hard. 1 || 3 || 4; or 1 || 3, 4 soft; or 1 ^3, 4 soft. 


> 


< 


< 


> 


1, 2, 4 hard. 1 || 2 || 4; or 2 || 4, 3 soft; or 2 ^4, 3 soft. 


< 


> 


< 


> 


All hard. 1 | 2, 3 | 4. Double-colUnear region. 


> 


< 


> 


< 


All hard. 2 | 3, 1 | 4. Double-colUnear region. 


< 


> 


> 


> 


All hard. 1 | 2. No contribution to NLO kernel. 


> 


< 


> 


> 


AH hard. 2 | 3. No contribution to NLO kernel. 


> 


> 


< 


> 


All hard. 3 | 4. No contribution to NLO kernel. 


> 


> 


> 


< 


All hard. 1 | 4. No contribution to NLO kernel. 


> 


> 


> 


> 


All hard. No singularities. No contribution to NLO kernel. 



Table 1: Classification of the four parton phase space. The notation < (>) means that the correspond- 
ing invariant is smaller (greater) than S min - 
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3 

X 52 H^g - ZZi)&{si2 - SMn)@{s23 - Srmn)&{su " Smin) 
(=1 

+ ^N^l ^//?'^(fci,^2,^3)c(°)(3,l,2)|A('')(l,2,3)p 

x[5(xg — xi) + 5{xg —X2) + 8(xg —X3)] 

X 0(521 - *min)0(*31 - >ymin)0(>5'23 " Smin) , (5. 19) 

where we have used the cyclic invariance of the color-ordered amplitudes together with the freedom to 
relabel the momenta. We adjusted the boundaries of the z-integrals as described in section |4] The last 
term will contribute only to 5- type terms; dropping it and including the contribution from 2 || 3, 3 || 4, 
4 II 1 we obtain, 

a, 1 



2N' [ dz[ dR'{l,2,3) |a(o)(1,2,3)|2^-/^(o)( 
Jo J 271 z 



X ^ 8(;Cg - ZZi)®{sn - Sram)®{s23 - Sirmi)®{sn - ^min) • (5.20) 
(=1 



We can further rewrite it using eqn. (I4.27I) . 



2N' j dz I dR'iW) |A(o)(l,2,3)|2£5(^,-x,-z)g-^(")(z) 
2N'J dz J dR\W) |A(0)(l,2,3)|2£6(^,-x,-z)g V)(z) 

-N^ I dR\\,2) [5(3;-xiz) + 6(y-X2z)]|A(°)(l,2)|22(5(°)(l,2) + c('')(l,2,l))^ V)(z) 
J 271 £ 

l/z^') ® - ( 0^(0) 0^(0) 

271/ £ V27I/ £2 

-N^JdR^iia) [5(y-xiz) + 5(3;-X2z)]|AW(l,2)|22(5('')(l,2) + c('')(l,2,l))gV)(z). 

(5.21) 

The last term does not show a proper factorization because the combination 5 ( 1 , 2) + C ^''^ ( 1 , 2, 1 ) is 
still process dependent, but it will eventually cancel against virtual corrections with similar structure. 

We continue with contributions where two invariants are smaller than ^min- There are two types, the 
third and fourth types in tabled The latter, 

&{Smm - ■5'12)0('5'23 " Sram)&{s^m " S34)&{si4 - Smin) 
+ &isi2-Smm)®{Smm-S23)&{s34-Srmn)®{Srr,m-Su), (5.22) 

correspond to configurations with two independent pairs of coUinear gluons. Here, the collinear inte- 
gral over the variable not present in a 8-function can be done, while the other remains, so we obtain 

JdzJ c//?f(l,2)c('')(z^i,^2,(l-z)^i)[5(xg-xiz) + 8(x,-x2z)]|A(i)(l,2)|2gl/:(o)(z), (5.23) 
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where we have adjusted the boundaries of the z-integral as usual, and have dropped terms which con- 
tribute only to the 5-function in the kernel. Note that the integral C {zk\ , 2, ( 1 — z)^i ) still depends on 
the momentum fraction z, 

c<».(*,fe.(.-*)-c<».(*„*„.,,^-Sf^(^)'(^)'i(.'+(i-.r-2). 

(5.24) 

In the third type of region in tabled two invariants with a common momentum ('adjacent' invariants) 
are smaller than Smin, 

©(■^min - Si2)®{s2-i - ^•min)0('5'34 " ■5'min ) © (■5'min " Su) 
+ ®{Smm - ^12)0(^min " S2-i)®{s?,'i " ■5'min)0('Sl4 " ■Smin) 
+ @{S12 - Smm)®{Sram " S2-i)®{Smm " ^34)0(^14 " ■^min) 

+ 0(^12 - 5min)0('5'23 " S,am)®{Smm " S34)@{Smm " Su). (5.25) 

These contain subregions with qualitatively different types of unresolved contributions, in which we 
must use different factorizations - as long as we are not using a limiting function unifying different 
limits (see for example ref. [69]). For example, consider the region defined by 

®{Smm - Su)®{Smm " ■5'23)0('S34 " Sram)®{su - ^min)- (5.26) 

Gluons 1, 3, 4 cannot be soft, otherwise the last two functions would vanish. The constraint may be 
satisfied in two distinct ways: gluon 2 can be soft, or three momenta can be collinear, 1 || 2 || 3 (with 
or without gluon being soft). 

To distinguish the different subregions, we may introduce additional functions. For example, multi- 
ply eqn. (I5.26t by 

1 =@{sm-Smm)+®{Smm-Sm){®{s2A-Smm)+®{Smm-S2A)) (5.27) 

with Sijk = {ki + kj + kkY to obtain 

©(■^■min - Sn)®{Smm " S23)@{sm " ■5'min)0('5'l4 " Sram)®{sm " ■S'min) 2 SOft 

+0(.?34 - ■Smin)0('Sl4 " ■Smin)0(*min " *24)0(*min " Sm) 1 || 3 and 2 SOft 

+0(^34 -5niin)0('5'l4-'5'min)0(*24-*min)0(*min-'5'l23)- 1 II 2 || 3 (5.28) 

In the first term, gluons 1,3, and 4 are resolved, while gluon 2 is soft, so this will contribute only to 
8-function terms. In the second term, the matrix element will not be singular enough to produce a pole 
unsuppressed by s„wi BtI . In the third term, we must use the triple-coUinear factorization. 

Before turning to the factorization and the computation of resulting integrals, let us consider the last 
type of region (the second type in table Q, where three nearest-neighbor invariants are smaller than 

'^min^ 

®{Smm - Sl2)®{Smm " S2-i)®{Smm " Sm)®{su " ■^min) 
+ 0('5'min - Sn)®{Smm " ■5'23)0('5'34 " ■5'min ) (■^'min " Su) 
+ 0('Smin - ■Sl2)0('S23 " ^mm)0(^mm " ^34)0(5min " Su) 

+ @{sn-Smm)®{Svmn-S23)®{Smm-S3A)®{Smm-SiA)- (5.29) 
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To match the term we selected above, examine 



0(^34 - ^min)0(^min " ■Sl4)0('ymin " Sn)&{Smm " S23)- (5.30) 

Once again we must introduce additional ©-functions to distinguish between configurations in which 
different factorization formulae apply. In the case at hand, we must distinguish between 

1,2 soft; lsoft,2||3; or 2 soft, 1 || 4. (5.31) 



Multiplying eqn. (I5.30t by 

1 = {&{S13 - Smin) +&{s^in- Sl3)){&{Synm - S24) +&{S24 - Smin)) (5.32) 

we obtain 

&{S34- - ■5'mm)0('5'min " ■5'23 ) © (■5'min " Su)&{Smm " Su) 
X ^0(5min-'y24)0('5min-'Sl3) 1, 2 SOft 

+ e{s24 - ^min)0(^mm " ^b) 1 SOft, 2 || 3 

+ 0(513 - SMn)@{Smm - S24) 2 SOft, 1 || 4 ) (5.33) 



where we have dropped the contribution containing 0(^13 — Smm)@{s24 — Smin) because it is kinemati- 
cally forbidden. The first term in eqn. (I5.33t contributes only to 5-function terms in the kernel, and we 
will not consider it further. 



We will want to attach the first of the remaining two terms to the region defined by eqn. (I5.28t . while 
the last term we will associate with a the similar product of theta functions with (1,2,3,4) —> (4, 1,2,3). 
To organize this, define functions aggregating various theta functions above, 

&cia,b,C;d) = &{Sad - Smm)&{sbd - SMn)&{Scd - Srmn)&{Sram- Sabc), 

&w{a,b,c) = &{Smm-Sab)®{Smm-Shc)&iSmin-Sac), 
&s{a;b,C;d) = &iSmm-Sad)&{sbd-Srmn)&iScd-Smm)- (5.34) 

The complete surviving contribution of terms with two or three small invariants, other than the double- 
collinear already accounted for in eqn. (I5.22t . is then given by inserting 

£ [0c(pi,p2,p3;p4) 

p=(1234), (2341), (3412), (4123) 

+ 0s(Pi;P2,P3;P4)0w(Pl,P2,P3) 

+ 05(p3;P2,Pi;p4)0iv(pi,p2,p3)] (5.35) 

into the integrand of the last term in eqn. (I5.lt . 

Each of the different permutations p corresponds roughly to a triple-collinear region, with added re- 
gions where one of the 'outer' gluons becomes soft. That is, they correspond to regions where a 
three-particle invariant becomes small. Let us focus on one of these contributions, say as above the 
region where 5123 0. 
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Using the factorizations I47ll48ll49l 

|a('')(1,2,3,4)|2 |SpIit"-'='=(l,2,3)p x |a('')(1 +2 + 3,4)|2 (5.36) 

and 

|a(°)(1,2,3,4)|2 ^'^"^M(lsoft,2 || 3) x +2 + 3,4)|2 (5.37) 

the integrand in the region where sns — > is given by 
Im = I (^1+^2 + ^3, 4) 1 2 

X ^05(l;2,3;4)0vi'(l,2,3)M(lsoft,2 || 3) 
+ 05(3;2,l;4)0w(l,2,3)M(3soft,l || 2) 

+ 0c(l,2,3;4)|Split°''='=(l,2,3)p^ + non-singular. (5.38) 

As noted in ref. BTll the function describing the mixed 'soft-collinear' Umit, M(asoft,& || c), can be 
derived from the triple-collinear splitting function. Consequently, the difference of the two limiting 
functions - the soft-collinear one and the triple coUinear one - gives only a finite contribution when 
integrated over a region where both are valid. In particular, as far as extracting poles is concerned, 
contributions of the form 

05(l;2,3;4)(M(lsoft,2 || 3) - | Split"'"'=(l,2,3)p) (5.39) 

can be dropped. We can use this to combine different contributions so as to simplify the structure of 
the resulting integrals we must compute. In general, the fewer constraints (theta functions), the better. 

Using 

©(■^34 - ^min)0(^14 " ^min)0(^24 " ^min) = 

1 - &{Smm - Su)&{s24 - ^•min)0('5'34 " Sram) " " ■5'min)0('5'24 " ■5'min)0(*min " ■5'34) 

-&{Smm - S14)&{S24 - Smm)@iSmm " S34) " " S24) (5.40) 

along with eqn. ( I5.39t . we obtain 

/l23 = (0(^min-^123)|Split'^'='=(l,2,3)P 

+ 05(l;2,3;4)(0„,(l,2,3)-0(5n,in-5i23))M(lsoft,2||3) 
+ 05(3;2,l;4)(0w(l,2,3) -0(5n,in -^i23))M(3soft, 1 || 2) 

- (^min -Su)&{s24- ^min) © (^min " ^34 ) © (Smin " ^ 123 ) | Spht"'='= ( 1 , 2 , 3 ) ^ 
-0(^min-^24)0(^min-^123)||Split'^'='=(l,2,3)|2^|A(")(l+2 + 3,4)|2. (5.41) 

The theta functions in the penultimate term require both gluons 1 and 3 to be soft, and hence that term 
can contribute only to 6-function terms in the kernel. The theta function in the last term forces gluon 2 
to be soft, and hence as discussed above, the term will not contribute any singular terms. 

One may also be tempted to drop the second and third terms. This temptation should be resisted, 
because although the region is small, the soft-collinear factorization function is nonetheless sufficiently 
singular to produce a contribution, as shown by a careful analysis in a different context fTOl . 
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Let us first evaluate the primary contribution, given by the first term in eqn. (I5.41t . The other regions 
related by cyclic invariance (where respectively S234, ^'134, and ^'124 vanish) give equal contributions; 
adding all four, we obtain 



24^4 f dR\l,2,3,4) |A(°)((l+2 + 3),4)|2£5(^,-xO|Split''--(l,2,3)|20(.^i„-.i23). (5-42) 



3 

L 

i=l 

Factorizing the phase-space measure in the triple collinear region fTOll . 



dR%\,2,3A) ^ —dR''{P,4)dRU{l,2,3), (5.43) 



with 



fif/?fon.(l>2,3) = r(l-2£) ^ ^dzidz2dZ3dsi23dsi2dsi3ds23 

5(1 - zi - Z2 - 23)8(5123 - si2 - si3 -S23), (5.44) 

A= {Z3S12 - Z1S23 - Z2Sl3f - 4ziZ2S23Sl3, (5.45) 

and 

ki=Zi{h+k2 + k3)=ZiP (5.46) 



the contribution of eqn. (15.421) takes the form, 

IN' J dR'{P,4)dRUX\,2,3) |a('')((1+2 + 3),4)|2 
3 

£ 8(^, - xi) I Split"-'='=(l , 2, 3) |20(^^i„ - ^123) 
i=l 



N'J dz J dR'{ia)\A^'\ia)\H^{^s-ziz) + 5{x,-z2z))[^yiK^'\z) 
^y'-[h^'^^K^\x,) (5.47) 



where 



2 1 3 

-Kj:'\z) =N^J J<n.(l,2,3)^Xj5(z-Z/)|Split''''='=(l,2,3)|20(.^i„-.i23). (5.48) 

Note that the factorization of the phase space given in eqn. (15.441) is exact up to an additional factor 
(1 — if one defines the 'momentum fractions' outside the collinear region by 



2ki ■ k/t 2k j -kA , ,9 

Zi= ~ ~ =^^, Sijti = {ki + kj + kk + kif (5.49) 

2kijk ■ K4 '^1234 



with j = 5i23Ai234 and 



h = -^—k^, kijk = ki + kj + kk- yh- (5.50) 
\—y 
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Using eqn. (I5.44t together with l47l . 

£ |Split'-'=(l,2,3)|2 = 2x| 

ph. pol. ^ 

2 



(1 -e) (22^123 - (1 -Z3)^23)" 2(l-£>23 , 3(1 -S) 



■^12'*123 (1-23)^ ^12^^23 ^4 



123 



^ 1 / (1-Z3(1-Z3))' ^ zi+Z2Z3+zi ^ Z2Zl-ziz3-2 ^ Z2 



■5'12'5'123 V Z3Zl(l-Zl) 1-Z3 Z3(l-Z3) 1 - Z3 

+ 1 /^3„2 3 (2-Zi+zf)(zi+Zi(l-Zi)) ^ 1 ^ 1 



2i'12^23 V Z3(l-Z3) Z3Z1 ( 1 - Z3) ( 1 - Zl ) , 

+ (■5'12 ^■5'23,Zl ^23), (5.51) 

(there is a factor of 1/4 included here compared to ref. f^T\ to account for our normalization of the 
color matrices) we obtain 

1 (-102^3 + 55z^ + 105z^ - 102z + 55) \ 1 ,on ^ ^^w ^2 

-;?(-z)52(z) - 8(1 +z)Li2(z) - 8/7(2) ln(z) ln(l -z) 

^2 , 1 73z^+z + 88 , 

-107?(z)ln(l -z)^ + ln(z) 

z 

1 -102z3 + 55^4 + 105^2 -I02z + 55, , 67,, 
+ 3 In(l-z)-yMz) 

1 -330z + 95 + 101z2 , 2 / / N . 4, 



36 (1-z) V ' 3 

with p{z.) defined in eqn. (14.141 . and where fl^ 



+ 71^ 7^(z) + -(l+z) , (5.52) 



^2(2) 



f i+z dw , 


1 — w 


/ —In 




w 


w 



-2Li2(-z) - 2 lnzln( 1 + z) + ^ ln2 z - J . (5.53) 
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In the derivation of eqn. ( 15.521 ) we have used the integrals given in the appendix. 

Next, we compute the two additional contributions from the second and third terms in eqn. ( 15 .4 U . The 
limiting function which we need to integrate is given by l47l 

M(/soft,i II k) = 2 ^'-'J + ^^"(W + ZjSu. + Sij)sj, + s„ ^ ^^^^^ 

with / being the momenta of the adjacent color connected hard parton in the antenna containing the 
soft gluon. (The normalization is again different from 1.47 J on account of different normalization con- 
ventions for the color matrices.) The term 

0s(l;2,3;4)(0w(l,2,3)-0(^™„-5i23))M(lsoft,2||3) 

= @{S34 - Srnm)@{s24 - Smm)&{Smm - Sl4){&w{^,2,3) - @{Sij,in - Si23))M{lS0ft,2 \\ 3) (5.55) 
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in eqn. (I5.41I) yields the following contribution. 



6gtN' f ^/?^(1,2,3,4)|a(0)(1,2,3,4)|2£6(x,-x,) 

1=1 

(■^34 - ■^min ) © (■^24 " ^min ) © (^min - 5l4)(0*(l,2,3)-0 {Smin - 5 1 23 ) ) M ( I SOf t , 2 1 1 3 ) 
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(5.56) 



where 



Ksiz) = UStz^N^ J t//?^„ii.(l,2,3)(6(z-Z2) + 5(z-(l-Z2)))M(lsoft,2|| 3) 

X 0(^34 -■5rain)0('524 -■5min)0('5min " ^m) (©W ( 1 , 2, 3) - ©(^min " ■5123))- (5.57) 

The easiest way to obtain ^5(2) is to calculate the contributions from 0vi?(l,2,3) and 0(smin — .S123) 
separately and then take the difference. The integration over the region given by ©(^'mm — '^123) is 
similar to that which yields eqn. (I5.48t . For the integration over the region given by 01^(1,2,3) we can 
use the result given in L70J . Subtracting the two contributions the final result reads: 



2e 



1 



SpA 

Sr^n) r(l-2£) 



Piz) 
iz) 



(5.58) 



£3 Vr(i-£)2^ r(i-£)2 

One can also compute the integral directly, as a check, and we obtain the same result. Including all 
other terms with the same structure (3 soft, 1 || 2 and cyclic permutations), we finally obtain 

N^J dz J dR'^{P,A) \A^^\l,l)\\h{xg-zxi) + h{xg-zx2)\Ki{z). (5.59) 



Combining the separate contributions from eqn. (I5.21t . eqn. (I5.23t . eqn. (I5.47t . and eqn. (I5.59t our 
final result for the singular part of the four parton final state is. 



g)' ^^(0) + l;^(o) ^ _ V) ^^(0) j ^ 



(5.60) 



in which, 



-2(5(°)(l,2) + c(°)(l,2,l))|V)(z) + c(^zl,2,(l-z)l)|l^(^z)+(|y^^^^ 



a,\2 

271/ 



271/ £2 



ni-£) 



A 7 \ 2e 



2izil-z)rpiz) 



-2-^£+(^7r2-^)£^ + 0(£3) 



(5.61) 



5.3. Assembling the Kernel 

Upon combining the results for the singular contribution of three- and four-parton final states, eqn. i5Jl 



and eqn. (15.601) respectively, we obtain 

^^(0) +;^(0) ^ ( 1 _ 1^(0) ^ l^iO) + 1 ^ ^ 

E \e e e e-^ e J 



(5.62) 
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Note that the equation above contains only and 1/s singularities, the poles present in indi- 
vidual terms cancel in the sum. In addition, the singularity is independent of as it ought to be. 
Inserting this result in eqn. ( I3.15I ) we obtain 



dxg 



e \e e e e-^ e J 

\ 8 8 2 8 8 4e^ 2 8 y 

8 

+ /j(o) ^ (^-^S,'P^'^ - 2^7^^) +0(8"). (5.63) 

Thanks to the leading order result in eqn. (I4.32t . we see that 

8 

is finite, so that we may write. 



dxg 



\8 8 28^ 8^48^ 2 8 / 

+ ®-(5eP(°)+7^W) 
8 

_ /j(0) ^ _L +2SeP('') ®/^(°)) +0(8°). (5.65) 

The term involving h^'^^^ is finite as it should be. The last term can be written in a more suggestive 
form, 

_/j(0) ^(5eP(°) +7^(°)) ® (SeP(°) +7^(°)). (5.66) 
The leading-order result implies that 

5eP('')+^(°' (5.67) 

is of order 8, and the convolution will not produce additional singularities. This has been checked 
via an explicit calculation. We may therefore drop the term in eqn. ( I5.66t . The requirement that 
the remaining, first, term on the right-hand side of eqn. (15.651) be finite then allows us to extract the 
next-to-leading order kernel P^^\ 

\s^P^'^ = (^-^(^)-i^(V^^(<')®^(")-^V(z) + ^5e2po/'(°') . (5.68) 

Intuitively, the first two terms are in some sense the radiative corrections to the leading-order spUtting 
amplitude, while the following two terms, 

l^(o)^^(0)_l^(^) (5.69) 

28 8 
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just remove the iteration of the leading-order kernel. The last term can be thought of as an ultraviolet 

(1) ^.(1) 



renormalization. Plugging in the explicit results for Kr , Kv , 1^ [eqn. (I5.15t . eqn. (I5.17t . eqn. (I5.52t . 



eqn. (I5.61t l. together with 
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441 



44 



= {16n^9{,Nf {5{\-zH'' + ^p{z) + 12-^-- Uz+^z 

V J 5 Z 

fl - 47 + 372 + 74) 
-81n(z)p(z) +4^ ]Z \ Hz) + 81n(l 



+ e 



2(l-£) 

9 z 



16(l+z)Li2(z) 



11 +2z+ llz^)ln(l -z) - 12;7(z)ln(l -z) 



3 z 

4 - ( 1 1 + 3z + 12z2 ) In (z) + 8 ( 1 + z) In (z)^ 
3 Z 

, ^0 11 11 / s 2t 67,,, 

-2(21n(z)' - y ln(z) - y ln(l - z) - -Ti' + -)piz) 



+ 0{e^), (5.70) 



we finally obtain, 



,./27,, , 67,2 1. ,11 25 44 1,,,, ,,,,2 

V2 9 z333z 

-,22 1 T 67 \ 

+ (41n(z)ln(l-z)-31n(z)2 + yln(z)--7i2 + y)/,(z) + 2p(-z)52(z)j (5.71) 



in agreement with known results for the time-like kernel finilT^IOI . 



6. Conclusion 

Intuitively, the Altarelli-Parisi kernel summarizes that part of collinearly unresolved radiation from a 
short-distance process which must be absorbed into the scaling evolution of descriptions of initial- or 
final-state hadrons. The computation described above makes this direct connection precise, and shows 
how to use it in order to compute the kernel. The approach described in the present paper effectively 
breaks down the NLO computation into smaller and simpler parts, whose intermediate terms have an 
independent meaning and are subject to consistency checks on their own. 

As described in section 15.31 the approach effectively computes the the next-to-leading order kernel as 
a radiative correction to the leading-order kernel (after cancellation of soft singularities), less an itera- 
tion of the leading-order kernel. As usual, the radiative corrections arise from a virtual conection and a 
real-emission correction. Both are computed from higher-order splitting amplitudes, quantities which 
govern the collinear behavior of gauge-theory amplitudes. These intermediate quantities are useful 
elsewhere in their own right, and can be checked independently. For example, the virtual corrections 
are basically given by the one-loop splitting amplitude, which should satisfy certain supersymmetry 
identities 154J . The real-emission correction arises from integrating the splitting amplitude describ- 
ing the behavior of tree-level amplitudes in a gauge theory as three color-connected partons become 
collinear. 
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The complete computation of the NNLO corrections to the Altarelli-Parisi kernel remains an important 
goal for particle theorists. These corrections, and parton distribution functions relying on them, are 
needed for programs evaluating jet production to NNLO at hadron colliders. In the approach described 
in this paper, the required ingredients would be the two-loop 1 — > 2 splitting amplitudes; the one- 
loop 1 — > 3 splitting amplitudes; and the 1 — > 4 splitting amplitudes. The second ingredient has been 
calculated by Catani, de Florian and Rodrigo ITTll and the third by Del Duca, Frizzo and Maltoni f72l . 
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thank Zvi Bern, Sven Moch, and Stefan Weinzierl for a careful reading of the manuscript. 



A. Double unresolved phase space integration 

The calculation of contributions from the four-parton final state in particular eqn. ( I5.48t leads to the 
consideration of integrals of the following form 



C„,.=.,r.f / j,..j,.,__(_A)-i- (A.1) 

with 



1 



-A>0 -' 



A = {sik{l - Zi - Zk) - ZiSjk - SijZkf - ^ZkZiSjkSij. (A.2) 
The integration can be done in <i = 4 — 28 dimensions yielding 

Ci,i = -2iilsTl-^hi%%'-''^{\-zd\l-z.kf2F,{-z,-z,l-z ), 

1 -l-2e„-E„-E„-E 1 



1 



co,i = -7I-^o•^'vz7V(l-Z'•)"^ 



and 



Z^C2,o -2zfo(l -Z.)C2,-l + (1 -Zc)'C2,-2 

- TT-^-^.-i-^^z-z-z-f 1 -z?(e- ^^^^ + ^^) + iii^±^^ (A 4) 

- \l-2e'«^^ Za) [e Zt+ZaZ^^izs+ZaZc?)- ^^'^^ 

Note that the C2.m integrals can not appear in arbitrary combinations in a gauge theory amplitude. The 
individual integrals appearing in eqn. iAAi are not regularized in dimensional regularization, in con- 
trast the specific combination is regularized. The combination given above is exactly the combination 
appearing in eqn. (I5.48t . 
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